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THE  TRANSIENT  RESPONSE  OF 
GELATIN  TARGETS  TO  PROJECTILE  IMPACTS 

ABSTRACT 

This  report  presents  results  obtained  in  the  course  of  an  investiga¬ 
tion  of  the  transient  response  of  gelatin  targets  to  projectile  impacts 
at  1.37  kilometers  per  second.  The  report  has  two  main  sections.  The 
first  describes  the  experimental  technique  used  to  gather  transient 
response  data  for  the  targets,  and  presents  and  discusses  the  test  results. 
The  second  section  is  devoted  to  development  of  a  data  reduction  method 
through  which  observed  particle  velocity-time  response  for  the  elastic 
deformation  of  a  target  may  be  analyzed  to  yield  load-time  response.  Some 
of  the  present  data  are  reduced  using  this  method  and  the  indication  is 
that  a  useful  supplement  to  the  observed  results  has  been  achieved. 
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I.  INTRODUCTION 


For  several  years  an  experimental  program  has  been  conducted  in  the 
Ballistic  Research  Laboratories  (BEL)  to  study  the  transient  response  of 
semi-infinite  targets  to  projectile  impact.  The  intent  of  the  program  is 
to  provide  data  that  will  assist  in  the  formulation  and  evaluation  of 
theoretical  analyses  of  the  hypervelocity  impact  problem1*.  Accordingly, 
an  induction-wire  technique  has  been  developed  to  record  the  velocity¬ 
time  response  of  particles  situated  on  the  (deformation)  axis- of -symmetry 
of  normally  impacted  target b2’  1+**.  The  technique  has  been  used  expen¬ 

sively  with  wax  blocks.  Collected  data  were  reduced  to  supply  the  speed 
of  propagation  of  the  shock  front  generated  in  the  targets  by  impacting 
pellets,  as  well  as  the  peak  particle  velocity,  normal  stress,  and  density 
ratio  at  the  front.  To  date,  however,  no  use  has  been  made  of  the  observed 
velocity  histories  subsequent  to  initial  particle  disturbance  by  the  shock. 
Significant  advantages  would  be  gained  if  methods  were  available  to  reduce 
these  data  in  terms  of  stress  levels,  strains,  and  deformation  rates. 

Toward  this  end,  a  reduction  technique  is  developed  in  this  report  for  the 
case  of  elastic  target  deformations.  Attention  was  restricted  to  the 
elastic  case  to  provide  a  basis  for  study  of  more  complex  situations 
involving  permanent  deformations. 

The  reduction  technique  to  be  presented  is  applicable  to  a  diversity 
of  problems  involving  high  intensity,  short  duration  loads  with  axial  or- 
spherical  symmetry.  It  is  most  clearly  given  in  terms  of  specific  data. 
Hence,  the  development  here  will  be  with  reference  to  experimental  results 
obtained  in  an  investigation  for  the  Hugoniot  of  a  gelatinous  target  mat¬ 
erial.  At  the  impact  velocity  used  (l.37  kilometers  per  second)  the 
targets  displayed  a  large  region  of  elastic  behavior;  therefore,  the  data 
are  convenient  for  present  purposes. 


*  - - 

Superscript  numbers  denote  references  found  on  page  37. 

-** 

The  technique  is  applicable  only  to  electrically  nonconducting  targets. 
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The  report  is  arranged  in  two  sections.  The  first  describes  the 
experimental  method  and  presents  the  data  along  with  its  reduction  to 
provide  the  characteristic  features  of  the  front  of  the  target  disturbance. 
This  includes  the  Hugoniot  of  the  material.  The  second  section  develops 
and  applies  the  method  for  analyzing  the  elastic  response  of  the  targets. 

II.  THE  EXPERIMENT 

1.  Description  of  Experiment  and  Presentation  of  Results 

The  experiments  consisted  of  firing  0.50-caliber  cylindrical  projec¬ 
tiles  3/8  inches  in  length  at  9-inch-cube  target  blocks.  The  gun  was 
smooth  bore.  The  experimental  data  collected  were  velocity-time  records 
for  particles  on  the  axis-of- symmetry  of  the  deformation  field  generated 
by  the  normal  impact  of  the  projectile  on  the  target  face.  The  projectile 
material  was  ethocel  (projectile  mass,  approximately  0.5  gfam)  while  the 
target  blocks  were  gelatin  (see  Appendix  A  for  the  casting  and  instrumen¬ 
tation  procedures  of  this  particular  material).  The  impact  velocity  for 
all  data  given  here  was  nominally  1.37  kilometers  per  second.  The  particle 

velocity-time  data  were  obtained  by  use  of  the  imbedded  wire  technique 

* 

developed  initially  for  use  in  studies  in  wax  .  Adaptation  for  use  with 
the  present  material  presented  no  practical  difficulties. 

In  order  to  have  a  sufficiently  complete  description  of  the  axial 
response  of  the  targets,  and  to  get  the  necessary  data  to  determine  the 
Hugoniot  curve,  the  wire  particle- velocity  monitors  were  placed  at  six 
stations  along  the  axis-of -symmetry  of  the  targets.  A  spacing  of 
1.27  centimeters  between  wires  was  used,  with  the  first  wire  located 
1.27  centimeters  from  the  target  face  (Figure  l).  A  representative  set  of 
oscilloscope  records  from  a  firing  is  shown  in  Figure  2.  These  records 
provide  a  velocity  history  of  the  particles  at  each  station,  and  time-of- 
arrival  data  for  the  front  initiated  by  the  impact.  With  the  arrival  data, 
the  average  velocity  of  propagation  between  wires  is  given  by 
_ 

Details  of  this  experimental  technique  have  been  presented  and  discus¬ 
sed  previously.  The  reader  is  referred  to  References  2,  3,  and  4  for 
details. 
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A  SECTION  THROUGH  THE  AXIS  OF  SYMMETRY  OF 
THE  TARGET  SHOWING  THE  LOCATIONS  OF  THE  PARTICLES 
MONITORED  BY  THE  WIRE  PICK-UPS 


TIME,  5  y.  SEC/CM  — « 


A  REPRESENTATIVE  SET  OF  DATA  RECORDS 
FOR  AN  ETHOCEL-GELATIN  IMPACT. 
IMPACT  VELOCITY  1.37  KM/SEC. 
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Fig  2 


(II-l) 


where  D  is  the  wire  spacing,  millimeters,  and  T  is  the  travel  time  between 
wires,  microseconds.  For  the  ethocel-gelatin  impacts  studied  here,  the 
velocity  of  propagation  data  are  shown  in  Figure  3  as  a  plot  against  dis¬ 
tance  from  the  point  of  impact.  Data  are  shown  only  far  travel  times 
taken  from  the  records  of  individual  oscilloscopes,  i.e.,  for  1.27  to 
2.5k  centimeters;  3.8l  to  5.08  centimeters;  and  6.35  to  7.62  centimeters. 
Although  all  three  oscilloscopes  are  simultaneously  triggered  at  impact, 
zero-time  reference  from  photograph-to-photograph  is  not  sufficiently  well 
determined  to  provide  accurate  elapsed  times. 

The  behavior  demonstrated  by  the  data  of  Figure  3  is  qualitatively 
as  expected.  The  front  of  the  disturbance  propagates  with  decreasing 
velocity  as  it  penetrates  the  target  block,  asymptotically  reaching  a 
value  of  1.55  kilometers  per  second.  This  occurs  by  the  time  the  front 
arrives  at  the  3.8l-centimeter  station.  The  upward  trend  of  the  curve 
between  5  and  7  centimeters  was  a  repeatable  characteristic  of  these  par¬ 
ticular  tests  and  was  found  not  to  be  an  instrumentation  error.  However, 
results  of  tests  not  reported  here  show  this  trend  is  not  to  be  expected 
to  continue.  It  is  most  probably  due  to  a  weak  nonlinear  characteristic 
of  the  target  material.  The  velocity  of  1.55  kilometers  per  second  may 
be  taken  to  be  the  elastic  dilatational  wave  velocity,  C^. 

Figure  4  presents  the  maximum  particle  velocity  in  the  deformation 
field  just  behind  the  advancing  front  of  the  disturbance.  These  data  are 
gotten  by  noting  the  peak  electromotive  force  on  the  oscilloscope  records 
of  the  type  shown  in  Figure  2.  The  voltage  is  directly  proportional  to 
the  particle  velocity^.  The  scatter  in  results  is  due  to  variation  in 
impact  velocity  and  other  round-to-round  irregularities.  It  is  more 
evident  for  the  particle  velocity,  u,  than  for  the  velocity  of  propagation, 
U,  shown  in  Figure  3  since  the  former  is  a  more  sensitive  function  of 
conditions  at  impact . 
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VELOCITY  OF  PROPAGATION  OF  THE  DISTURBANCE 
FRONT  AS  A  FUNCTION  OF  DISTANCE  FROM  THE  IMPACT  POINT 
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MAXIMUM  PARTICLE  VELOCITY  VS  DISTANCE  FROM  IMPACT 
FOR  ETHOCEL-GELATIN  AT  1.37  KM/SEC 


The  velocity  of  propagation,  U,  and.  maximum  particle  velocity  behind 
the  disturbance  front,  u,  given  in  Figures  3  and  ^4- ,  respectively,  are 
sufficient  information  to  determine  the  Hugoniot  curve  for  the  target 
material.  The  data  points  for  this  curve  are  obtained  through  use  of  the 
conservation  relations 

aR*  =  -pq  Uu  (II-2) 


T 

T 

O 


(H-3) 


where  a*  =  maximum  radial  stress  at  the  front  of  the  disturbance 

R 

p  =  undisturbed  density  of  the  target  material,  1.02  grams  per 
°  cubic  centimeters 

t  pn 

—  =  —  =  specific  volume  ratio  at  the  front  of  the  disturbance. 

To  p 

Figure  5  is  a  plot  of  the  Hugoniot  curve  obtained  from  the  data  shown  in 
Figures  3  and  4.  Up  to  a  stress  level  of  4  X  10®  dynes  per  square  centi¬ 
meter,  the  data  are  well  represented  by  a  straight  line  whose  slope  gives 
a  value  for  U  of  1.55  kilometers  per  second.  Above  this  load  level,  the 
curve  shows  a  slight  but  definite  increase  in  slope.  This  was  to  be 
expected  from  the  behavior  of  U  at  points  closer  to  the  target  face  than 
3.81  centimeters. 

To  extend  the  Hugoniot  curve  to  stress  levels  greater  than  achieved 
here  will  require  further  experiments.  It  is  worth  noting,  however,  that 
the  target’s  ma.ior  constituent  is  water.  Gross  differences  between  the 
Hugoniots  for  the  two  materials  are  not  to  be  expected.  Within  the  load 
range  covered  here,  this  is  evidenced  by  noting  that  the  propagation 
velocity,  U,  for  water  is  about  1.45  kilometers  per  second. 

The  terminal  appearance  of  the  targets  provides  a  useful  supplement 
to  the  transient  data.  A  hemispherical  region  of  intense  cracking  and 
fracture  exists  around  the  impact  point.  The  radius  of  this  region  is 
2.5  to  3.0  centimeters.  A  few  isolated  fracture  surfaces  usually  extend 
to  depths  of  5.0  centimeters;  however,  there  is  no  other  evidence  of 


COMPRESSIVE  STRESS, -Or,  DYNES/CM  xIO 


permanent  deformation  beyond  3.0  centimeters  from  the  target  fa.ce.  Hence, 
the  indication  is  that  the  targets  behave  elastically  past  a  depth  of 
roughly  3-5  c  ent imet  er  s . 

2.  Discussion 

A  brief  description  of  the  response  of  the  targets  is  now  possible. 
This  will  be  useful  in  the  analysis  of  the  particle  velocity  records.  The 
pertinent  results  are  given  by  Figures  2,  3,  4,  and  5,  and  the  terminal 
target  appearance.  The  impacting  cylinder  produces  (and  sustains  for 
short  period)  an  intense  load  pulse  that  propagates  into  the  target.  As 
this  proceeds,  its  amplitude  attenuates  with  a  corresponding  decrease  in 
the  velocity  of  propagation.  Eventually,  at  a  depth  of  about  3-5  centi¬ 
meters,  the  pulse  is  reduced  to  the  level  of  elastic  deformations  and  the 
signal  thereafter  propagates  at  the  elastic  dilatational  wave  velocity. 

The  kinetic  and  potential  energy  deposited  in  the  material  near  the  impact 
point  produces  continued  deformation  and  fracture  of  the  target  for  a 
considerable  period  after  passage  of  the  shock  front.  No  details  of  this 
process  are  provided  by  the  present  data. 

The  specific  characteristics  of  the  particle  velocity-time  behavior 
are  shown  by  the  records  in  Figure  2.  From  the  impact  point,  only  the 
positions  at  3.81  centimeters  and  greater  are  of  interest  as  this  is  where 
elastic  deformations  occur.  Response  of  the  particles  is  qualitatively 
the  same  at  each  station.  The  material  is  undisturbed  until  arrival  of 
the  front  which  is  signaled  by  a  rapid  rise  in  the  particle  velocity.  A 
peak  is  attained  in  approximately  6  microseconds.  Afterwards,  in  about 
20  microseconds,  the  velocity  decays  to  an  appearant  asymptotic  value. 
Actually,  this  is  the  beginning  of  a  very  gentle  decrease  to  zero  velocity 
that  lasts  several  hundred  microseconds. 

Station-to-station  attenuation  in  pulse  intensity  is  evidenced  by 
the  decrease  in  peak  particle  velocity  plotted  in  Figure  4.  A  similar 
attenuation  exists  in  the  asymptotic  value  reached  by  the  particle  veloc¬ 
ity  after  the  peak.  The  latter  decreases  with  distance  from  impact 
approximately  as  the  inverse  square  of  distance. 
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The  rapid,  rise  in  particle  velocity  with  the  arrival  of  the  disturb¬ 
ance  front  indicates  the  front  has  the  nature  of  a  shock  wave.  However, 
the  resolution  of  measurement  is  adequate  to  show  distinct  structure 
within  the  shock;  that  is,  a  definite  rise  time  (6  microseconds)  to 
maximum  velocity  is  evident,  as  is  a  strongly  damped  high  frequency 
(250  kilocycles)  oscillation  during  the  initial  rise.  This  oscillation 
produces,  in  most  cases,  a  stepped  appearance  to  the  response.  Through 
examination  of  the  pickup  and  monitoring  circuitry  shows  this  oscilla.ting 
behavior  is  not  electronic  in  origin.  It  may  be  reasonably  presumed  to 
be  a  characteristic  of  the  target  response  to  the  conditions  imposed  by 
impact . 


III.  ANALYSIS  OF  VELOCITY  RECORDS 
1.  Preliminary  Assumptions 

The  aim  of  this  analysis  is  to  use  the  available  particle  velocity¬ 
time  records  to  deduce  further  characteristics  of  the  target's  elastic 
behavior,  e.g.,  stresses,  deformation  rates,  etc.  Two  assumptions  are 
necessary  to  obtain  this  information.  The  first  is  that  at  the  positions 
and  times  of  interest  the  response  is  linearly  elastic.  The  stipulation 
of  linearity  presumes  the  material  obeys  Hooke's  law  in  the  elastic  region, 
that  deformations  are  small,  and  that  time  effects  on  the  elastic  behav¬ 
ior  are  negligible.  The  magnitude  of  the  deformations  may  be  checked  a 
posteriori;  however,  other  deviations  from  classic  elastic  behavior  are 
not  so  easily  ascertained  without  extensive  material  testing.  This  does 
not  seem  necessary  as  it  may  be  noted  that  the  behavior  demonstrated  by 
Figure  3  does  not  indicate  any  strong  nonlinear  effects. 

The  second  assumption  is  that  along  the  axis  of  symmetry  the  target 

response  depends  only  weakly  on  the  coordinate  6  (see  Figure  6)  and, 

therefore,  may  be  considered  spherically  symmetric.  This  assumption,  in 

conjunction  with  the  one  of  elastic  behavior,  has  been  used  by  other 

investigators0  in  the  study  of  axially  symmetric  problems  with  high- 

intensity,  short  duration  loads.  Data  were  interpreted  by  comparison  with 

analytical  solutions  to  specific  problems.  In  contrast,  the  intent  here 

is  to  employ  the  assumptions  to  develop  a  direct  method  of  data  reduction. 

IT 
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WAVE  STRUCTURE  FOR  AN  IMPULSIVE  LOAD  APPLIED 
NORMAL  TO  THE  SURFACE  OF  A  SEMI-INFINITE  SOLID 


19 


ASSUMED  WAVE  STRUCTURE 


The  assumption  of  spherical  symmetry  results  in  considerable  simpli¬ 
fications  in  the  analysis  of  data;  however,  certain  of  its  limitations 
must  be  kept  in  mind.  An  immediate  implication  is  that  no  rotational 
(shear)  disturbances  are  present  in  the  deformation  field.  However,  a 
general  feature  of  elastic  disturbances  involving  axial  symmetry  and  con¬ 
centrated  loads  is  the  presence  of  both  shear  and  head  waves  following  a 
forerunning  dilata.tional  wave.  An  analysis  neglecting  their  presence  can 
be  used  only  so  long  as  the  deformations  are  reasonably  assumed  free  of 
rotational  disturbance.  The  time  interval  for  which  this  is  the  case  may 
be  estimated  by  reference  to  Figure  6a,  (a  meridian  plane  representation 

of  the  wave  pattern  associated  with  an  impulsive  load  F  normally  applied 

8 

to  the  surface  of  a  semi-infinite  elastic  solid  ).  Shown  are  the  fronts 
of  the  dilatational,  shear,  and  head  waves  initiated  by  application  of 
the  load  at  time  t  =  0.  The  disturbance  in  region  ' 


is  dilatational 


while  that  in  "b"  is  a  combined  dilatational  and  shear.  A  typical  parti¬ 
cle,  Q,  on  the  axis-of- symmetry  OZ  is  undisturbed  until  arrival  of  the 
P-wave.  Between  this  arrival  and  arrival  of  the  S-wave,  the  deformation 
at  Q  is  irrotational.  The  elapsed  time  between  P  and  S  arrivals  and  the 
time  for  which  the  spherically  symmetric  representation  might  be  valid, 
therefore,  is 


T  =  0Q 


('s "  0  ’ 


(III-1) 


where  Cg  is  the  velocity  of  shear  waves,  and  is  the  velocity  of  dila¬ 
tational  waves.  Applying  this  representation  to  the  present  case,  the 
formula  given  by  Equation  (ill-l)  serves  to  estimate  the  time  interval 
for  which  an  analysis  based  on  spherical  symmetry  may  be  applied.  For  the 
gelatin,  Cp  =  1.55  kilometers  per  second  while  Cg  may  be  taken  to  be 
roughly  Cc  =  0.3  kilometer  per  second*;  therefore,  with  0Q  =  3.8  centi- 

O 

meters,  Equation  (ill-l)  gives 

i  , 


T  «  3.8(33  -  6)  =  102  psec. 


*  This  results  from  an  assumed  value  for  Poisson*s  ratio,  v,  of  0.^8. 
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This  indicates  the  spherically  symmetric  assumption  may  he  used  without 
concern  for  effects  of  the  shear  wave  for  about  100  microseconds  at  the 
3. 8l- centimeter  station  and  proportionally  longer  for  the  succeeding  sta¬ 
tions  .  In  the  following,  no  attempt  will  be  made  to  analyze  data  records 
for  more  than  30  microseconds  after  arrival  of  the  dilatational  wave. 

A  further  comment  with  regard  to  shear -wave  effects  must  be  made. 

In  Figure  6b,  a  representation  more  nearly  depicting  the  real  situation 
than  in  Figure  6a  is  shown.  Elastic  behavior  holds  outside  the  shaded 
region  defined  by  R  =  R^j  inside,  combined  effects  of  strong  nonlinearity, 
plastic  behavior,  and  bifurcation  are  anticipated.  The  response  within 
R  =  R-^  defines  the  boundary  conditions  for  the  elastic  region.  It  is  to 
be  noted  that  the  previous  considerations  include  the  implicit  assumption 
that  these  boundary  conditions  generate  a  P  and  a  S-wave  structure  simi¬ 
lar  to  that  shown  in  Figure  6a.  No  direct  verification  of  this  is 
available. 

It  is  recalled  that  the  assumption  of  spherical  symmetry  was  taken 
on  the  basis  of  weak  dependence  on  8  of  the  target  response  along  its 
axis  of  symmetry.  Some  justification  of  this  is  desirable  even  though  it 
is  plausible.  Also,  inquiry  should  be  made  into  the  adequacy  of  a  spheri¬ 
cally  symmetric  representation  of  the  weakly  8 -dependent  situation. 

These  are  considered  in  turn. 

Direct  experimental  or  analytical  evidence  is  not  available  to 

justify  the  assumption  of  weak  9 -dependence.  However,  recourse  is  made 

to  investigations  of  two  related  problems,  the  results  of  which  indicate 

7 

the  assumed  behavior.  The  first  is  that  of  Rae  and  Kirchner  and  concerns 
the  initial  stages  of  the  process  of  hypervelocity  impact.  For  this 
problem,  they  have  concluded  9 -dependence  of  the  deformation  field  on  the 
symmetry  axis  is  weak  and  that  a  spherically  symmetric  analysis  adequately 
represents  the  axial  behavior.  The  second  result  is  the  elastic  solution 
available  for  a  point  normal,  impulsive  force  applied  to  the  surface  of  a 
_ 

T  is  quite  large  owing  to  the  gelatin  having  a  Poisson*  ratio  of  nearly 
0.5.  For  more  common  values,  the  interval  T  will  be  considerably 
reduced. 
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semi-infinite  solid.  Here,  it  is  found  that  the  6 -distribution  of  ampli¬ 
tude  coefficients  for  the  P-vave  (see  Reference  8)  is  weakly  0-dependent. 
These  two  problems  bracket  the  current  one  in  the  severity  of  the  loading 
conditions  imposed  at  the  body  surface.  Their  demonstration  of  the  weak 
0 -dependence  lends  substance  to  assuming  it  here. 

Final  consideration  is  given  to  the  adequacy  of  the  spherically 
symmetric  solution  to  represent  the  weakly  0-dependent  problem.  Interest 
is  confined  to  distances  from  the  impact  point  that  are  only  several  times 
greater  than  R^,  the  boundary  of  the  damaged  region.  Figures  6b  and  6c 
represent  the  problem.  The  latter  figure  gives  a  comparison  of  0- 
independent  and  0 -dependent  conditions  on  R  =  R^,  At  0  =  0,  R  =  R^,  they 
are  identical  and  therefore  both  present  acceptable  solutions.  The  ques¬ 
tion  is  whether  this  continues  for  R>  R^(9  =  0).  The  affirmative  is  seen 
to  hold,  provided  R  and  the  time  t  are  not  too  large,  since  the  response 
at  a  typical  point  Q  is  an  integral  of  the  effects  produced  by  the  loading 
over  the  entire  surface  R  =  R^.  The  effect  produced  by  a  load  applied  at 
0  =  «/ 2  is  felt  at  Q  at  a  time  later  (in  proportion  to  the  distances 
A Q,  BQ)  than  one  simultaneously  applied  at  0  =0.  Furthermore,  the  mag¬ 
nitude  of  the  response  produced  by  loading  at  B  is  attenuated  more 
severely  than  a  corresponding  signal  from  A  due  to  BQ  being  greater  than 
AQ. 

2.  Analysis 

It  is  well  known  that  elastic  disturbances  possessing  spherical 
symmetry  may  be  represented  by  a  single  potential  function  satisfying  the 
spherical-wave  equation,  that  is, 

fc(R,t)  =  ^  '  (HI-2) 

.  ,  R 

where  T  =  t  -  — 

Cp  =  velocity  of  elastic  dilata.tional  waves 
t  =  time 

R  =  distance  from  the  center  of  symmetry. 
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The  single  displacement  component,  uR,  is  given  by 


UR  =  "ciR  ’ 


(III-3) 


and  the  stress  components,  crR,  are  provided  by  Hooke’s  lav  as 


E 


aR  (1  +  v)(l  -  2vj 


du. 


a  -  v)  ^  ri  - 


u. 


R 


E 


°e  (1  +  v)(i  -  2v) 


UR  ^R| 

T  +  v  ‘STl  ' 


(m-4) 


(III-5) 


vhere  E  and  v  are  Young’s  modulus  and  Poisson’s  ratio,  respectively.  The 
corresponding  strains  are 

(m-6) 


^R 
CR  ~  "§T” 


e6  R  ‘ 


(III-T) 


Within  the  framework  of  Equations  (III- 2)  through  (ill  7),  it  is 
seen  that  the  assigned  task  of  determining  characteristics  of  the  defor¬ 
mation  or  load  fields,  in  addition  to  the  observed  velocity  history  at  a 
specific  axial  station,  reduces  to  determination  of  the  function  F(t) 
[Equation  (111-2)3 .  Two  approaches  to  this  problem  are  to  be  developed 
here.  The  first  involves  empirical  construction  of  a  function,  F(t),  that 
adequately  represents  the  velocity  history  at  a  specific  particle  position 
r  =  R  .  This  function  is  then  used  to  determine  the  stress  and  velocity 
behavior  within  the  elastic  region  of  the  target.  The  second  approach 
involves  a  graphical  analysis  of  the  individual  records  at  each  station, 
r  =  R*,  to  provide  the  function  F(t  -  B*/^)  for  these  stations  separately. 
With  this  method,  it  is  possible  to  evaluate  the  desired  response  charac¬ 
teristics  at  each  station  from  the  particle  velocity  history  at  that 

station. 
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The  first  of  the  two  approaches  is  appealing  owing  to  the  analytical 
compactness  with  which  it  represents  the  target  deformation.  Although 
guidance  is  provided  by  known  solutions  to  dynamic  problems  having  spheri¬ 
cal  symmetry,  a  major  drawback  to  the  method  is  the  effort  required  to 
establish  the  "adequate"  representation  of  the  available  data  record. 

The  second  approach  does  not  possess  this  difficulty  and  has  the  attribute 
of  being  direct  and  devoid  of  ambiguity.  It  is  the  more  appropriate  of 
the  two  for  reduction  of  data  from  several  stations,  as  in  the  present 
ca.se,  or  for  interpreting  large  amounts  of  data..  The  two  methods  are  con¬ 
sidered  separately  below. 

a .  Construction  of  a  Function  F(t)  to  Represent  an  Observed  Velocity- 
Time  Record.  The  characteristic  behavior  of  the  velocity-time  records 
has  been  discussed  previously.  For  purposes  of  analysis  attention  is 
restricted  here  to  interpretation  of  the  records  from  a.  single  test,  the 
round  shown  in  Figure  2.  Also,  selection  of  the  function  F  is  based  on 
the  data  obtained  at  a  single  station,  R  =  5-08  centimeters.  The  adequacy 
of  the  selected  representation  may  then  be  checked  by  determining  whether 
it  satisfactorily  portrays  the  velocity  history  at  the  other  data  stations. 

Enumeration  of  all  the  details  associated  with  selecting  the  specific 
function  used  to  represent  the  velocity  record  would  be  difficult.  Briefly, 
it  was  noted  that  the  impulsive  nature  of  the  projectile  loading  was  in 
many  ways  similar  to  a  uniform  load  suddenly  applied  to  the  interior  of  a 
spherical  cavity.  The  response  to  such  a  loading  situation  consists  of 
exponentially  damped  circular  functions.  With  this  guide,  a  satisfactory 
function  F  was  found  to  be  of  the  form 

F  =  Fl(t)  +  F2(t)  ,  (m-8) 


2? 


/  -1  k 

where  F^t)  =  e  sin(£  t  +  a)  , 


F2(t)  =  -A2  {(sin  a2)t  + 


+  |)2  cos(|32  t  +  a2)]j 


F1  =  F2  =  0  for  r  <  0  , 


cot  a,  =  -r—  . 
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The  function  describes  the  high  frequency  oscillations  while  F2 
describes  the  gross  behavior  of  the  deformation  response.  From  Equa¬ 
tion  (III-2),  it  is  to  be  noted  that  t  =  0  represents  arrival  of  the 
disturbance  at  any  station,  R. 

Applying  Equation  (III-8)  to  Equations  (III-2)  and  (III-3),  the 
velocity  history  is 


B  B  +  F2 

V  =  V1  +  V2  =  'St  - R - 


1  1 


■V  r  [K2  -  K  M£l1/2 


cosC^  T  +  a±  +  p±) 


i[(7i2  +  ^12)]  sin  h  T} 


(III-9) 


v2  sin  o:2  x  -72t  f  y^  /^\2 

A2  ~  R2  "RS  VR"Cd/+\Cd/ 


sin(t>2  t  +  a2  +  P2) 
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where 


v  -  h 
tan  P1  =  ■  afffi  ' 

Vcd 

tan  02  "  (l/R  -  72/Cd)  • 

A  restriction  placed  on  Equation  (III-9)  is  that  v  he  zero  for  T  =  0. 

This  lea.ds  to  the  relation 

4(,12  -  t)x2)sin  ax  -  27 ±  ^  cos  aj 

+  A sin  a2  -  P2  cos  a£j  =  0  .  (HI-10) 

The  parameters  in  Equation  (III-9)  that  define  a  specific  profile  of 
velocity  versus  time  are  7  ,  7 2,  and  The  7's  axe  decay  factors 

while  the  f)*s  establish  frequency  of  oscillation.  The  constants  and 
A2  are  amplitude  coefficients  for  the  two  functions  and  vg.  Appropri¬ 
ate  values  of  the  profile  factors  for  the  R  =  5.08-centimeter  curve  of 
Figure  2  are: 

7^  =  0.800  microsecond  1  (psec 
^  «=  2.520  psec-1 
72  =  0.145  psec  1 
P2  =  0.069  psec  1. 

Due  to  the  rapid  attenuation  of  v^,  the  amplitude  coefficient  Ag  is 
established  with  sufficient  accuracy  by  requiring  (v2)/max  at  R  = 

5.08  centimeters  be  equal  to  the  maximum  velocity  from  the  data  record. 
This  results  in 

_  7  7 

Ag  =  64.5  X  10  in  /sec 

which, when  applied  to  Equation  (III-IO),  gives 

A1  =  1.01  X  10"10m5  . 
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Figure  7  is  a  plot  of  the  velocity  history  at  R  =  5-08  centimeters 
as  transcribed  from  the  data  record  and  as  given  by  Equation  (III-9)  with 
the  above  stated  values  of  7^,  y  ,  etc.  A  good  description  of  the  actual 
record  is  seen  to  have  been  obtained  except  for  details  of  the  initial 
oscillating  response.  A  better  fit  could  have  been  obtained  by  further 
trial  and  error  adjustment  of  the  parameters;  however,  the  gain  achieved 
would  not  warrant  the  effort. 

Now  that  the  specific  function  F  has  been  established  to  represent 
the  velocity  behavior  at  R  =  5*08  centimeters,  it  may  be  used  to  evaluate 
other  response  characteristics  in  the  target.  First,  a  comparison  should 
be  made  of  the  predicted  and  observed  velocity  response  at  the  remaining 
data  stations.  This  is  also  shown  in  Figure  7.  The  representation  of 
actual  data  is  seen  to  be  excellent  at  all  the  stations  for  the  times 
(t  ^  0)  of  interest.  The  agreement  provides  some  support  to  the  assump¬ 
tion  of  spherical  symmetry.  If  the  assumption  were  grossly  unreasonable, 
marked  discrepancies  between  the  observed  velocity  response  and  that 
provided  by  F  should  result. 

Next,  the  function  F  is  used  to  examine  further  characteristics  of 
the  target  response.  Most  significant  of  these  are  the  normal  stresses, 
ctr,  Oq,  and  the  maximum  shear  stress  ( crR  -  aQ ) / 2.  These  are  obtained 
directly  by  application  of  Equation  (lII-8)  to  Equations  (III-4)  and 
(lll-5)j  with  the  appropriate  values  of  y^,  y ....  This  has  been 
done  for  the  radial  stress  cm,  and  maximum  shear  stress  (a,.,  -  aQ)/ 2.  The 
results  are  shown  in  Figures  8  and  9- 

The  stress  cu,  has  a  profile  in  time  very  similar  to  that  of  the 
x\ 

velocity  at  each  station.  It  quickly  reaches  a  peak  (details  of  the 
oscillatory  behavior  have  been  omitted)  at  about  6  microseconds  and  then 
decays  to  what  appears  to  be  a  residual  value  of  the  order  of  30  percent 
of  the  maximum  stress.  This  shows  the  highly  transient  character  of  the 
large  radial  stresses  imposed  by  the  impact.  The  present  data  indicate 
the  duration  of  intense  loading  to  be  about  15  microseconds.  A  similar 
circumstance  may  be  assumed  to  hold  in  the  region  of  damage,  R<  3.81  cm 
(centimeters) . 
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v  GIVEN  BY  F|  AND  F2 
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PARTICLE  VELOCITY  VS  RETARDED  TIME 
AT  THE  MONITORING  STATIONS 
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RADIAL  STRESS,  <TR,  VS  TIME  AS  GIVEN  BY  F2(Tj 


MAXIMUM  SHEAR  STRESS  VS  TIME  AS  GIVEN  BY  F2(t),  V0.48 


Figure  9  shows  the  behavior  of  (crR  -  Oq)/2,  the  maximum  shear  stress. 
(Again  the  initial  oscillatory  behavior  has  been  omitted.)  Here,  a  quali¬ 
tatively  different  behavior  is  seen  than  for  cr  .  The  shear  stress  reaches 
a.  peak  in  the  same  time  as  the  radial  stress,  but  no  significant  subse¬ 
quent  decay  occurs. 

b.  Graphical  Data  Reduction.  The  method  to  be  considered  now  will 
make  use  of  the  individual  velocity-time  records  at  each  station  to  deter¬ 
mine  additional  characteristics  of  the  load  process  at  these  stations. 

With  data  available  at  several  stations  as  in  the  present  case,  this 
provides  a  more  direct  means  of  data  reduction  than  the  previous  technique. 
Also,  it  may  be  used  when  only  a  single  data,  record  is  at  hand  to  provide 
information  for  the  region  near  the  data  station. 

Advantage  is  taken  again  of  the  assumption  of  spherical  symmetry. 

* 

Also,  use  is  made  of  the  conditions  uR=v=0atT=0  .  The  result  of 
these  conditions  is  that  F(o)  =  F;(o)  =  F  (o)  =  0. 

For  a  fixed  station  R  =  R*,  the  displacement  component  u^  is  readily 
obtained  by  a  first  quadrature  of  the  data  records;  that  is, 

ud(R*,t)  =  (  v(R*,£)d£,  t  >  0  (Ill-il) 

R 

while  u  (R*, t)  =0;  t  ^  0.  Applying  this  result  to  Equation  (ill- 5)  gives 
R 

FlLO  +  Ilii  =  -u(R*,t)  ,  (III-12) 

R  CD  R*^ 


where  primes  denote  derivatives  with  respect  to  the  argument,  T.  Multi¬ 
plying  Equation  (III-12)  by  the  integrating  factor  exp(CD/R*)T  gives 


d_ 

dr 


=  -(R*Ce) 


(C  /R*)t 
e  u(R*, 


T 


/ 


which,  upon  integration,  yields 


This  means  account  will  be  taken  of  the  structure  of  the  shock  wave. 
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F(t)  =  e 


(III-13) 


■(Cd/R*)t 


-(r*c. 


t  (C  /R)£ 

D  u(t)d| 


Equation  (III-13)  is  the  desired  result  .  It  allows  evaluation  of  F(t) 
at  any  of  the  stations  R  =  R*  by  a  simple  process  of  graphical  quadratures. 
Consequently,  any  desired  characteristic  of  the  target  response  is  imme¬ 
diately  available  through  a  direct  process  of  reduction  of  the  available 
data.  For  instance,  the  radial  stress  aR  is  given  as 


(i  +  u)(i  -  aj_l  (R*  T)  =  _(i  .  v) 

E  R 


R*‘ 


D  J 


-(Cd/R*)t  pT  (C-n/R*)! 
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'D' 
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.(£)a.6 


-  (l  -  *>) 


R* 


(III-14) 


Evaluation  of  the  necessary  integrals  in  Equation  (III-14),  that  is. 


uJR*,t)  =  v(R»,6)d6 

R  JQ 


and 


-  f 


■t  (cd/r*)i 


ur(r*,0^  , 


has  been  carried  out  for  the  data  shown  in  Figure  2  at  the  3.8l-,  5-08-, 
6.35-  and  7.52-centimeter  stations.  The  velocity-time  records  were 
transcribed  to  graph  paper,  and  a  planimeter  was  used  for  integration. 

The  results  obtained  are  shown  in  Figure  10.  Comparison  with  the  results 
provided  by  the  empirical  choice  of  F(t)  yields  agreement  within  10  per¬ 
cent.  The  graphical  method  provides  the  details  of  the  high  frequency 
motion  which  were  omitted  in  Figure  8. 


For  a  similar  development  involving  the  acoustic  wave  equation,  see 
Reference  3. 
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RETARDED  TIME,  T,  j-L  SEC 

RADIAL  STRESS,  Op,  AS  GIVEN  BY  GRAPHIC  EVALUATION  OF  F(t) 


The  most  attractive  feature  of  the  graphical  method  is  its  simplicity 
and  directness.  The  data  at  each  station  can  he  independently  reduced. 
Also,  it  is  to  he  noted  that  this  technique  provides  a  very  useful  pro¬ 
cedure  in  tests  where  velocity-time  data  are  available  from  only  a  single 
record.  The  functions  F(t),  F7(t)  can  he  obtained  graphically  from  rela¬ 
tions  such  as  Equations  (III-15)  and  (III-12)  at  the  data  station.  These 
functions  may  then  he  used  for  other  values  of  R  provided  proper  account 
is  taken  of  the  retarded  time 

R  -  R* 

T  =  t - p -  , 

4) 

and  as  long  as  the  assumption  of  spherical  symmetry  is  valid. 

IV.  SUMMARY  AND  CONCLUSIONS 

1.  Determination  of  the  Hugoniot  curve  for  the  gelatin  target 

2 

material  has  been  accomplished  for  stress  levels  up  to  10  dynes  per 
square  centimeter.  In  this  load  region,  the  curve  differs  little,  hut 
distinctly,  from  that  for  water.  For  the  load  intensities  covered  by  the 
present  tests,  this  result  was  to  he  anticipated  .  At  higher  load 
levels  the  same  situation  may  he  expected.  Due  to  the  practical  interest 
for  use  of  this  particular  material,  a  definitive  determination  of  the 
Hugoniot  is  desirable  for  a  more  extensive  range  of  loads  than  covered 
here.  Tests  are  in  preparation  for  such  an  extension  and  will  he  pre¬ 
sented  subsequently. 

2.  The  curves  of  Figures  2,  3,  and  h  summarize  the  characteristics 
of  the  front  of  the  disturbance  traveling  into  the  target.  The  data  pro¬ 
vided  furnish  a  detailed  account  of  initial  features  of  the  signature 
imposed  on  the  block  by  the  projectile.  A  discussion  of  these  features 
was  given  in  Section  II-B  for  the  present  impact  situation.  The  ability 
to  examine  these  features  is  significant.  For  many  problems  of  practical 
interest  involving  impacts  at  velocities  attained  by  service  weapons,  the 
evaluation  of  effects  of  projectile  geometry  is  important.  Appreciable 
under  standing  of  these  effects  could  be  provided  through  systematic  study 
of  the  initial,  transient  characteristics  they  impose  on  targets. 
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5.  The  data  analysis  presented  in  Section  III  of  the  particle 
velocity-time  records  represents  a  significant  extension  of  the  basic  data 
provided  by  the  inductance  wire  pickup  technique.  Considerable  attention 
has  been  given  to  justification  of  the  necessary  assumptions  for  making 
this  extension.  This  practice  seems  to  have  been  neglected  by  other 
investigators  in  the  past.  The  information  provided  by  the  data  analysis 
affords  an  opportunity  for  more  detailed  experimental  investigation  of 
impact  phenomena  than  has  been  possible  in  the  past.  Immediate  applica¬ 
tion  to  the  study  of  hypervelocity  impact  phenomena  is  possible  as  would 
be  its  use  in  further  studies  in  gelatin.  It  should  also  be  pointed  out 
that  the  basis  of  the  techniques  developed  here  is  equally  applicable  to 
various  problems  involving  high  intensity  loading.  The  most  notable 
problems  are  geophysical  and  involve  phenomena  such  as  earthquake  behavior 
and  surface  or  underground  explosion  effects. 

The  foundations  and  application  of  the  data  reduction  technique  must 
be  considered  preliminary.  Establishment  of  its  accuracy  and  applica¬ 
bility  can  result  only  from  further  study  and  use. 
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APPENDIX 


Preparation  of  Gelatin  Targets 

The  target  material  used  for  these  studies  was  a  mixture  of  collagen 
and  water.  Generally,  lots  of  thirty  pounds  were  prepared  in  a  ratio  of 
four  parts  water  to  one  part  collagen  by  weight.  The  preparation  proce¬ 
dure  was  as  follows: 

1.  Collagen  added  to  water  at  l80°F . 

2.  Constituents  mixed  for  3  hours  at  l80°F. 

3.  Liquid  gelatin  cast  into  9-inch  cubes  and  allowed  to  gel  at 
room  temperature  overnight. 

k.  Cast  targets  chilled  to  40°F  in  refrigerator  for  24  hours. 

5.  Targets  removed  from  refrigerator  and  tested  with  a  maximum 
expo sure  of  20  minutes  to  room  temperature  of  68  F. 

Method  of  Wire  Placement 

The  inductance-wire-particle-velocity  pickups  were  emplaced  by  means 
of  a  rig  of  fine  thread.  A  pair  of  taut  piano  wires  was  used  to  assist 
in  alignment  of  the  0.003-inch  copper  monitor  wires.  These  were  removed 
prior  to  gelling  of  the  blocks.  A  plane  view  of  the  arrangement  is 
depicted  in  Figure  11. 


Swift's  USP  Grade  Pharmagel  A. 
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r  UPPER  SURFACE  OF  GELATIN 


A  CROSS  SECTION  OF  THE  TARGET  MOLD 
SHOWING  THE  METHOD  OF  WIRE  PLACEMENT 
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